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We review a number of ideas related to area law scaling of the geometric entropy from the point of view 
of condensed matter, quantum field theory and quantum information. An explicit computation in arbitrary 
dimensions of the geometric entropy of the ground state of a discretized scalar free field theory shows the 
expected area law result. In this case, area law scaling is a manifestation of a deeper reordering of the vacuum 
produced by majorization relations. Furthermore, the explicit control on all the eigenvalues of the reduced 
density matrix allows for a verification of entropy loss along the renormalization group trajectory driven by the 
mass term. A further result of our computation shows that single-copy entanglement also obeys area law scaling, 
majorization relations and decreases along renormalization group flows. 



PACS numbers: 03.75.Ss, 03.75.Lm, 03.75.Kk 

I. INTRODUCTION 

The amount of entanglement present in a quantum state is 
of fundamental relevance to determine how hard it is to simu- 
late it by classical means. It is generally argued that a highly 
entangled quantum state carries a huge superposition of prod- 
uct states that cannot be handled on a classical computer. Yet, 
this statement must be made precise, since a small amount of 
entanglement can indeed be simulated efficiently. The rele- 
vant precise question is, thus, how much entanglement can be 
efficiently simulated classically. 

This abstract question should at least be clarified when con- 
sidering relevant physical systems. Can the amount of en- 
tanglement present in a two-dimensional lattice of harmonic 
oscillators be efficiently represented in a classical computer? 
Although the answer to this question is not yet settled, quali- 
tative progress has been recently achieved. One of the ingre- 
dients essential to this discussion is the area law for the ge- 
ometric entropy and the representation of quantum states by 
projected entangled pairs. 

A important related problem is to understand how entan- 
glement varies along renormalization group (RG) trajectories. 
We shall bring growing evidence for the idea that RG flows 
entail a loss of entanglement. This entanglement loss will be 
shown compatible with area law scaling of the entropy. 

We organize the contents of this paper by first reviewing 
a number of previous results on area law scaling of ground 
state entropy in different systems using the language of con- 
densed matter, quantum field theory and quantum information 
theory. We shall then present a computation of entanglement 
entropy on a discretized bosonic free field theory in arbitrary 
dimensions. This gives us control on the eigenvalues of the re- 
duced density matrix on a subsystem which, in turn, allows for 
a discussion of majorization relations obeyed by the reduced 
density matrix of the system. We extend this discussion to the 
single-copy entanglement measure. RG loss of entanglement 
is also verified in detail for arbitrary dimension networks of 
harmonic oscillators. 



II. A BRIEF REVIEW OF THE AREA LAW 

A. Measures of entanglement for many-body systems 

An arbitrary quantum state of e. g. N spins is, in general, 
highly entangled. To quantify such a statement we can use 
various figures of merit. For instance, concurrence Q| is eas Y 
to compute and detects some pairwise entanglement though it 
cannot scan correlations throughout the system. An appropri- 
ate and widely used candidate to quantify entanglement is the 
von Neumann entropy of the reduced density matrix of the 
state under analysis, when only a subset of degrees of free- 
dom is retained. To be precise, let us consider a quantum state 
made out of N qubits \tp) <G (C 2 )® N and its density matrix 
of a p = Next, we consider the reduced density ma- 

trix of a subset of qubits denoted by A, p A = Tr^p where 
all qubits but those belonging to the set A are traced out. The 
von Neumann entropy of the reduced density matrix is then 
defined as 

S(pa) = -Trp A logp A . (1) 

The entropy is often referred to as entanglement entropy. 

In general a set of particles will be distributed randomly 
over space. Entanglement entropy can be computed for all 
sorts of partitions of the system, yielding information about 
the quantum correlations among the chosen subparts. A par- 
ticular and extremely relevant class of physical systems are 
those made of local quantum degrees of freedom which are 
arranged in chains or, more generally, in networks. For such 
systems it is natural to analyze their entanglement by studying 
geometrical partitions, that is, computing the entanglement 
entropy between a set of contiguous qubits versus the rest of 
the system. We shall referred to this particular case of entan- 
glement entropy Qlas geometric entropy (also called fine 
grained entropy in 01). 

The appearance of scaling of the geometric entropy with the 
size of the sub-system under consideration has been shown to 
be related to quantum phase transitions in one-dimensional 
systems, further reflecting the universality class correspond- 
ing to the specific phase transition under consideration 0,0] 
(see also 10]). Broadly speaking, a large entropy is related to 
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the presence of long distance correlations, whereas a small en- 
tropy is expected in the presence of a finite correlation length. 
The precise scaling of geometric entropy does eventually de- 
termine the limits for today's efficient simulation of a physical 
quantum system on a classical computer. 

There are many other ways to quantify entanglement. The 
von Neumann entropy we have chosen as our central figure 
of merit has an asymptotic operational meaning. Given in- 
finitely many copies of a bipartite quantum state, it quantifies 
how many EPR pairs can be obtained using local operations 
and classical communication. A different measure of entan- 
glement can be associated with the analysis of entanglement 
on a single-copy of the quantum system. This single-copy en- 
tanglement ll7ll8ll9l floll can be defined as 

E x {p A ) = -log p$ (2) 

where p 1 ^ is the maximum eigenvalue of p A - This quantity 
provides the amount of maximal entanglement that can be ex- 
tracted from a single copy of a state by means of local opera- 
tions and classical communication. As we shall see later on, 
the von Neumann entropy and the single-copy entanglement 
appear to be deeply related in any number of dimensions. 



B. Volume vs. area law 

Random states are known to carry large entanglement. To 
be precise, let us consider a random infinite system o qubits. 
On average, the density matrix for a random subset of N 
qubits carries maximum von Neumann entropy, 

S(p N ) ~ N . (3) 

This result fllll shows that the entropy of random states grows 
as the number of particles included in the subset. This is re- 
ferred to as a volume law scaling. An arbitrary state uses the 
maximum possible superposition of the basis elements with 
no symmetry whatsoever among their coefficients. Its efficient 
representation by classical means appears certainly difficult. 

Physical theories create entanglement through interactions, 
which are typically local. Thus, e.g. the ground state of a 
sensible physical Hamiltonian is not a random state. It is nat- 
ural to expect a low amount of entropy since local interac- 
tions will entangle the non-contiguous degrees of freedom in 
a somewhat sequential way. We may encounter local intense 
entanglement that dilutes at long distance. This is precisely 
the structure of standard quantum theories, with correlations 
that decay with a power law at phase transitions and with an 
exponential law away from them. It is then reasonable to ask 
what is the limit of efficient simulability in terms of the entan- 
glement present in a given state. 

In many physical theories, local degrees of freedom are ar- 
ranged in a specific geometrical way as mentioned previously. 
We may have quantum systems defined on spin chains, net- 
works or, in general, D-dimensional lattices. Those systems 
may have a continuum limit described by a quantum field the- 
ory or, alternatively, may be devised as quantum simulators, 



a preview of quantum computers. We may then discuss the 
amount of geometrical entanglement present on the system 
from three complementary points of view: condensed matter, 
quantum field theory and quantum information. 

As we shall see, the basic ingredient of locality of interac- 
tions suggests that entropy for a geometrical region should be 
dominated by the entanglement present on the surface sepa- 
rating it from the rest of the system. To be precise, consider 
an infinite D-dimensional lattice where we assign part A to 
an inner hypercube of size L, N — L D , and part B to the 
outside. Locality seems to suggest 

S(p L ) ~ L°- x ~ N 2 ^ . (4) 

This behavior is commonly referred to as area Jaw scaling for 
the geometric entropy. Let us note that one-dimensional quan- 
tum systems correspond to a well understood limiting case for 
the above formula, where the power law turns out to be sub- 
stituted with a logarithmic scaling at phase transitions, that is 

S{p L )~\ogL, (5) 

and saturates away from them 

S(pl) < constant ,VX , (6) 

as shown in Ref. These results are deeply con- 

nected to conformal symmetry and control the classical simu- 
lability of the system. 

Recent evidence hints at a log violation of the area law 
in some two-dimensional systems made with anticommuting 
variables fl4lll5lll6lfl7ll . To be precise, some of these models 
display an entropy scaling law of the type 

S(p L ) -L*- 1 log L. (7) 

It is unclear whether such systems support a limiting quantum 
field theory description in the continuum limit. 

It is important to make a general remark concerning the 
different approaches to the computation of entanglement in 
quantum systems. Let us note that discretized quantum sys- 
tems allow for uncontroversial computations of the entropy. 
This is not the case of quantum field theories, where regular- 
ization and renormalization are needed since the number of 
degrees of freedom is formally unbounded. In such a frame- 
work, the adimensional entropy requires the appearance of 
some short-distance regulator e 

S(p L ) ~ \-j (8) 

which entails the necessary discussion of its renormalization 
and its observability. Let us just mention here that the coeffi- 
cient of the area law is universal for D = 1 systems whereas 
remains scheme-dependent in higher dimensions. 

The problem turns extremely subtle in the case of gravity, 
where the geometry of space-time is dynamical and the way 
to compute for a black hole the Bekenstein area law pre-factor 
from first principles is far from clear(H Qjl |2jJ. Recent 
progress on the side of AdS/CFT correspondence seems to 
link entanglement entropy in a quantum field theor y living on 
the boundary to the black-hole entropy of the bulk 119112011 . 
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C. Locality and PEPS 

The basic heuristic argument for an area law scaling of en- 
tropy for the ground state of physical systems is rooted in the 
locality of the interactions. Steps to make this ar gum ent quan- 
titative have been made in Refs. 12 ll l22l l23l l24l 1251 12611 . 

A local Hamiltonian tends to entangle nearest neighbors. 
Long-distance entanglement emerges as a coherent combina- 
tion of local interactions. The correctness of this argument 
would imply that the reduced entropy of a geometric bipar- 
tition of a system will get its main contribution from the en- 
tanglement between degrees of freedom at opposite sites of 
the boundary that separates the regions. This, in turn, implies 
an area-law scaling. Let us note that such a naive argument 
works in any dimension and does not depend on the corre- 
lation length present in the system. Area law would emerge 
from locality, whatever the mass-gap is. We shall discuss the 
limitations of this argument shortly. 

This argument needs a clear formulation and verification. 
Although we lack definite answers about the necessary and 
sufficient conditions a Hamiltonian must obey to produce 
a ground state with area law entropy, some progress has 
been achieved using one-dimensional Matrix Product States 
(MPS) and their generalization to higher dimensions, Pro- 
jected Entangled Pair States (PEPS). We first consider a one- 
dimensional system with open boundary conditions described 
by a MPS 

w = E ■ ■ ■ A t-i) • • • *»> . w 

where the sum extends to ii, . . . , i n = l,...,d, which 
are physical indices attached to local Hilbert spaces, and 
ai,..., q„_i = 1, ... x, which are ancillae indices. The ten- 
sors A l a p can be viewed as projectors from the ancillae indices 
to a physical one. This representation provides the basis for 
the density matrix renormalization group technique. 

The generalization of the MPS construction to higher di- 
mensional networks carries the name of PEPS. In a D- 
dimensional network, where ancillae degrees of freedom are 
linked to their nearest neighbors, the role of the MPS projector 
is taken by a tensor of the form 

Kip , do) 

where the physical indices span a £>-dimensional lattice and 
ancillae run from 1 to x- Again, the role of each tensor A is 
to project maximally entangled pairs connecting local neigh- 
bors onto a physical local space. Entanglement is thus carried 
by the links connecting ancillae. Each entangled pair, that 
is, each sum over one ancilla index hides a connecting bond 
of the type X)«=i ~Jx\ aa )- ^ one °^ ^ e two anc iU ae m the 
bond is traced out, the entropy for the remaining ancilla is 
S = logx- 

We are now in a position to present the argument in Ref. 
ll2lll showing that finite x PEPS entail area law scaling for 
the entropy. Let us assume that the ground state of a quantum 
system is described by a PEPS with finite x- It follows that the 



entropy of a subpart of the system is bounded by the number 
of bonds which are cut by the separating surface times the 
entropy per broken bond. This amounts to an area law 

S{pa) < (#cut bonds) logx ~ Arealogx. (11) 

A violation of the area law within the PEPS representation 
requires infinite-dimensional ancillae. 

We should again distinguish the one-dimensional case, 
where the ground state of infinite critical systems are known 
to carry logarithmic entropy @|, 

S( PL )~^log-, (12) 
6 a 

where a is the lattice spacing and c the central charge that 
characterizes the universality class of the phase transition. 
Yet, the boundary of a one-dimensional block is made by two 
single points. Such a state with logarithmic entropy cannot 
be represented using finite dimensional MPS and we must re- 
sort to arbitrarily large x- This limitation is at the heart of the 
problems that the DMRG technique encounters when applied 
to quantum phase transitions. On the other hand, the entropy 
is bounded away from critical points and MPS provide an ef- 
ficient way to represent the system. MPS states with finite x 
are often referred to as finitely correlated states. 

Coming back to higher dimensions, it is then a major is- 
sue to establish whether finite x PEPS can describe faithfully 
the ground state of physical systems. The fact that PEPS with 
finite x can incorporate an area law is appealing. Recently, 
a particular class of finite PEPS has been constructed that 
display polynomial decay laws, that is long range correlation 
I2al . These PEPS are also shown to describe ground states of 
frustration-free Hamiltonians and such states can approximate 
exponentially well any finitely correlated state. It is still un- 
clear whether the ground states of standard quantum systems 
fall into this description or, alternatively, they need infinite x- 
This may set apart what is efficiently simulable from what is 
not. 



D. Renormalization group transformations on MPS and PEPS 
and the support for an area law 

We have argued that one-dimensional finite x MPS can sup- 
port a maximum amount of entropy independent of the size of 
the system and that, in contradistinction, finite D-dimensional 
PEPS can accommodate an area law. Let us give an indepen- 
dent quantitative argument for this statement. 

Consider a renormalization group transformation of a MPS 
state with constant A defined by the coarse graining of two 
sites ll27ll 

min(d2,x2) 

A%A^A%= ]T hU^V^ (13) 
1=1 

where we have decomposed the product of two adjacent ma- 
trices using a singular value decomposition. We can under- 
stand the unitary matrix U as a change of basis on the new 
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coarsed degree of freedom and construct a new MPS with 
A'^ = A;V^L. Therefore, the ancillae indices close under 
such operation whereas the physical index grows. Upon it- 
eration of this operation, the range of the physical index will 
reach a maximum value \2 and will get locked to that value. 
This is the magic of one dimension. The long-distance proper- 
ties of the system are completely described by a single square 
effective matrix! Entropy is then bounded. 

The analogous argument in two-dimensional systems fol- 
lows a slightly different path. The coarse graining step reads 

A" A b i -ob 

"Z " »2>f> = A cd , (14) 

- Ei=i XlUi v « a > v ;,w (15) 

As before, we can absorb the global U as a change of the local 
coarse-grained basis and assign a new PEPS to Ay. Note the 
different growth of indices. On the one hand, physical indices 
merge in groups of four and would naively need a volume law 
increase, dA. On the other hand, the ancillae rank increase 
from %4 to (x4)2, that is, it follows an area law. Given that 
the singular value decomposition will be locked by the small- 
est dimension of the two above, the area law will define the 
rank of the tensor that contains the effective long distance de- 
scription of the model. The argument generalizes to D di- 
mensions where the PEPS A l a a D with a physical index 
i = 1, . . . , d and ancillae indices ai, . . . , ct2D = 1, ■ ■ ■ , X- 
A renormalization group transformation of this PEPS makes 
the new collective physical index to run %' = 1, . . . , d 2 , 
that is, with a volume law, and the new collective ancillae 
a'x , . . . , a' 2D — 1, • • • , x 2 > tnat is, as an area law. The sin- 
gular value decomposition makes all the long-distance prop- 
erties of the state to be contained in an effective PEPS with 
a number of degrees of freedom that grows with just an area 
law. The rank of the effective PEPS is log \ef f = 2 £) ~ 1 log x- 
From this simple argument, it follows that PEPS can support 
an area law scaling for the geometrical entropy. 



E. Some explicit examples of area law 

There is an extensive literature on computations of the en- 
tropy for particular cases that cannot be faithfully summarized 
here. 

One-dimensional spin systems (e.g. quantum Ising model, 
XX model and Heisenberg model) obey a logarithmic scal- 
ing at the critical point g Q2 Qj |2fJ |2§|3J . Away from 
the quantum phase transition point, the entropy gets saturated. 
This explicit computation falls into the universal scaling pre- 
dicted by conformal invariance. This result has been further 
verified and extended to many other quantum systems in one 
dimension. I 

The literature on computations of entanglement entropy in 
higher dimensional systems is far less extensive due to the 
difficulty to produce explicit results. The first analysis of the 
entanglement entropy in two- and three-dimensional systems 



Spin chains away from criticality S ~ constant 
Critical spin chains S ~ log N 

D-l 

D-dimensional harmonic networks S ~ N~o~ 
NP-complete problems S ~ N 



were done in discretized approaches to quantum field theory 
GIB EH- Further analysis showed that the entanglement en- 
tropy is related to the trace anomaly in curved space times 
giving an explicit relation between the actual results for free 
fermions and free bosons l32ll . 

Rigorous computations in discretized harmonic networks 
proved no departure from the area law 123113311 . Further anal- 
ysis of entanglement entropy on higher dimensional networks 
has been done in Refs. IMbI I^. 



F. Exceptions to the area law 

We have argued that the area law is deeply connected to 
locality of interactions. It is, therefore, reasonable to expect 
violations of such scaling in models with non-local interac- 
tions. This is the generic case of a quantum computation of an 
NP-complete problem. It has been numerically verified that 
this is the case when an adiabatic quantum computation is ap- 
plied to the NP-complete Exact Cover problem, a variant of 
the 3-SAT problem. Along the computation, the ground state 
becomes maximally entangled, that is its entropy scales as the 
volume of the system l37ll38ll . A physical quantum computer 
will definitely need to face the challenge of maintaining those 
huge fine-tuned superpositions of states. 

Locality of interactions is not the only ingredient that con- 
trols entropy. Entropy is related to the eigenvalues of the 
Schmidt decomposition of a system in two parts. If the sub- 
systems retain a lot of symmetry, the sub Hilbert spaces or- 
ganize themselves in representations of the symmetry group. 
This entails a reduction of the Schmidt number of the above 
decomposition, that is, a lower entropy. Such a counter mech- 
anism to reduce the entropy in highly connected systems has 
been explicitely checked in the case of the Lipkin-Meshkov 
model which is defined by a spin system fully and symmet- 
rically connected. Although it is tempting to argue that the 
system is infinite-dimensional (the geometry of the Hamilto- 
nian corresponds to a simplex of N — > oo vertices), the en- 
tropy scales only logarithmically, which is the actual bound 
for symmetric spaces l39ll . This logarithmic scaling of the en- 
tropy follows the one-dimensional log law, which might just 
be an accident. 

It should not come as a surprise that slightly entangled 
states that do not correspond to an eigenstate of a given Hamil- 
tonian dynamically evolve to highly entangled states under its 
action. This has been analyzed in Refs. 14(1 14111 even for sim- 
ple Hamiltonians like the quantum Ising chain. No area law is 
expected for slightly entangled random states when they are 
evolved with local Hamiltonians. 

As mentioned previously, a case of non-trivial violation of 
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the area law was first considered in HH and then analyzed in 
Some two-dimensional systems with anticommuting 
variables were found to display a log correction to the area 
law, that is, S ~ LlogL. On the other hand, some previ- 
ous computations for free Dirac fermions seem to produce no 
area-law violation in any number of dimensions. 

This issue deserves further investigation. Finally let us men- 
tion that in the computation of quantum corrections to the en- 
tropy of a black hole, logarithmic corrections have also been 
obtained l43ll . 



G. Physical and computational meaning of an area law 

We can attach physical meaning to an area law scaling of 
entropy in different but related ways. We may argue that en- 
tropy is a measure of surprise due to quantum correlations and 
that a state that obeys an area law carries less correlations than 
a random state. As the size of the inner block increases, we 
only get a reduced amount of surprise, compare to the maxi- 
mum possible, when discovering that our block was correlated 
to the exterior. It is then arguable that the theory that has pro- 
duced such a state may accept a simpler description. In some 
sense, this argument is implicit in the holographic description 
of some quantum systems. 

From a computational point of view, low entropy means 
small quantum correlations, that is, small entanglement. It 
is known that states that are only slightly entangled can be 
efficiently simulated by classical means |46]. A fundamental 
question is thus formulated: what entropy growth law can be 
efficiently simulated by a classical computer? 

So far, this question can only be answered partially. In 
one dimension, D = 1, quantum critical phenomena show 
a logarithmic scaling which cannot be reproduced using finite 
MPS techniques. Formally, the simulation remain efficient in 
the sense that to reproduce critical behavior we need \ to be 
polynomial in L. This, though, produces an obvious practical 
computational slowing down and limitation. A new promis- 
ing idea to represent a quantum state with a different and non- 
local tensor structure has been proposed in Ref. Wh with 
the name of multi-scale entanglement renormalization ansatz 
(MERA). The basic idea is to substitute a linear MPS repre- 
sentation with a RG-inspired construction that also identifies 
the key use of disentangling operations for blocks before pro- 
ceeding to a coarsed description. 

The question in two dimensions has been addressed in 1I2U1 
in a sequential way. A PEPS is taken as lines of spins that are 
collected into effective degrees of freedom which are further 
treated in a MPS manner. 



III. AREA LAW IN D DIMENSIONS 
A. The Hamiltonian of a scalar field in D dimensions 

Let us consider the theory of a set of harmonic oscillators in 
D dimensions which is expected to verify area law scaling of 
the entropy. A number of non-trivial issues can be discussed 
in this explicit example. First, we shall analyze the regular- 
ized version of a scalar free field theory in order to get its 
reduced density matrix when an inner geometrical ball is inte- 
grated out. Its eigenvalues can, then, be used to compute the 
geometrical entropy that will scale as dictated by the area law. 
Second, we can compare the behavior of the entropy to the 
one of the single-copy entanglement. Third, we can analyze 
whether area law scaling is backed by a deeper sense of order, 
namely majorization theory. 

Our computation will generalize the one presented in Ref. 
|01 to D dimensions. Let us consider the Klein-Gordon 
Hamiltonian 

H = \Jd D x (tt 2 (x) + + M 2 l<^)| 2 ) , (16) 

where ir{x) is the canonical momentum associated to the 
scalar field <fi(x) of mass /j. The Z?-dimensional Laplacian 
reads 

^ r D-i Q r dr^) r 2 ^ ' ■* 

where r = | a?| and L 2 is the total angular momentum operator 
in D dimensions. It is convenient to introduce the real spheri- 
cal harmonic functions Zu m \, which are eigenfunctions of L 2 
with eigenvalues 1(1 + D — 2). The set of numbers {m} stand 
for other Casimir and component labels in the group SO(D). 
We now project the angular part of the scalar fields tt and <fi, 

n{m}{r) = r~ J d D xZ i{m} (6i, ■ ■ ■ ,6 D _ 2 ,tp)7r(x) 

(18a) 

0i{™}( r ) = r ~ j dDxZ l{m}(Sl, ■ ■ ■ ,0d-2,<P)<I>(x) , 

(18b) 

where r, Q\ , ■ ■ ■ , 6d-2 and ip define the spherical coordinates 
in D dimensions. The Hamiltonian now reads 

l{m} 

where, 

An ultraviolet regularization of the radial coordinate in the 
above Hamiltonian will transform the scalar field theory into 
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a chain of coupled harmonic oscillators. This is achieved by 
discretizing the continuous radial coordinate r into a lattice of 
TV discrete points spaced by a distance a, 



1 N 
} = ^E( 7r Hm} 



i=i 

( ■ , \d-1 ( _ <fo{m},j 



1(1 + D -2) 



(21) 



The size of the system is L = (N + l)a, where a and L act as 
an ultraviolet and infrared cutoff respectively. We can com- 
pare this expression with the Hamiltonian of an open chain of 
N coupled harmonic oscillators, 



N N 



2 —*—V"3 



(22) 



and identify as 

1(1 + D -2) 



1 - 



r 
i 



+ p?\5, 



1 + 



2j 



D-l 



D-l 



i X i 



i + 



D-l 



D-l 



h3 



(23) 



where is the step function. 

B. Geometric entropy and single-copy entanglement 



We now proceed to trace out an inner geometric ball around 
the origin to obtain the reduced density matrix of the ground 
state of the system on the exterior of that ball. Following sim- 
ilar steps as in ||2|] we define 57 as the square root of K, that 
is K = 0,2. The gaussian ground state of the system can be 
expressed as, 



Mxi,-,x N ) =ir- N/4 (detn) 1/ ' i e 



(24) 



where x = (xi, xn)- We construct the density matrix p ou t 
by tracing over the inner n oscillators, 

p out (x, x') ~ e -h(* T T*+*"-v*')+* T -P-*' , (25) 



where (3 and 7 are defined by 



1 



/3=-B T A- 1 
2 



D 



(26a) 
(26b) 



and A = fl(l+n, l+n), B = n(l+n, n + 1+N) and C = 
fl(n + 1+N, n + 1+N) are sub-matrices of fL 

We proceed with the diagonalization of this structure ro- 

rp — 1/2 

tating and rescaling the variables x = V j D y where 
7 = V t ^idV and jd is diagonal. Using this transformation, 



7 becomes identity, — ► /?' = j D ' VPV 1 7 £) ±/ ^ and the 
density matrix reads 

Pout(y,y')~e-^ 2+ v' 2 ) + y T -?'-y' . (27) 

If we do the appropriate change of coordinates y = W ■ z 
(where W is an orthogonal matrix) such that W T ■ j3' ■ W 
becomes diagonal with eigenvalues (3[, we get p out as a ten- 
sor product of the two coupled harmonic oscillators density 
matrices, 



-1/2 



N—n 



Pout 



m~ n e 



i( z i+ z ?)+p[ z i z 'i 



(28) 



We can now compute the entropy associated to the reduced 
density matrix p ou t- This entropy can be expressed as a sum 
over contributions coming from each term in the reduced den- 
sity matrix tensor product structure, 



N-n 



~>l{rn} 



y] £>l{m},i(£i) 1 



(29) 



i=i 



where 



Sl{ m },i(^i) = -log(l - £l{ m },i) 



1 - 6{m}, 



■ log 6{m}, 



(30) 

is the entropy associated to each sub-density matrix in the 
product shown in Ea. (l28> and £ii m \j is the parameter that 
generates the eigenvalues of these densities matrices. Note 
that each eigenvalue £ = £j{ m -w entails a set of probabilities 
of the form 



Pn = (l-m n n = 0,l,2,3, 



(31) 



defined by & = $/(l + (1 - /3f ) 1/2 ) for each l{m} set. 

To compute the total entropy, we have to sum over all pos- 
sible values of {m} and I. 



S 



l{m} 



(32) 



We realize from Eq.(|3DJ that H^ m j only depends on I, so 
the entropy associated to its ground state will also be {m} 
independent, and therefore 



S = J2v(l,D)S l , 



(33) 



1=0 



being v(l, D) the degeneracy of the total angular momentum 
operator L 2 for a fixed I. In three dimensions, for example, 
{m} = m can go from — I to I so that v(l,Z) is 21 + 1, The 



7 



same computation in D dimensions requires the computaiton 
of the degeneracy of SO(D) representations 



u(l,D) 



l + D-l\ (l + D-3 
I )~{ 1-2 



(34) 



Given the explicit knowledge of all the eigenvalues of the 
reduced density matrix, we can also obtain a formula for the 
single copy entanglement Eq.(0. The largest eigenvalue of 
density matrix for two coupled harmonic oscillators is (1 — 
£). This largest eigenvalue of the density matrix p ou t will be 
the product of the largest eigenvalues of the density matrices 
which compound p ou t, 

N-n 

pint = n n ^ ~ ^{m},*) 

l{m} i=l 
oo /N-n 

= n n (! ^w.o (35) 

;=0 \ i=l / 
The single copy entanglement finally reads 

oo /N-n \ 



(36) 



1 = 



i=i 



C. Perturbative computation for large angular momenta 

Note that our expressions for the entropy and the single 
copy entanglement depend on a final sum that ranges over all 
the values of angular momentum This sum may not be con- 
vergent as the radial discretization we have implemented is 
not a complete regularization of the field theory. To be pre- 
cise, the asymptotic dependence on I should be under control 
in order to correctly assess the convergence of the series. 

Let us note that, for I 3> N, the non diagonal elements of 
K Eq.j23> are much smaller than the diagonal ones. These 
suggests the possibility of setting up a perturbative computa- 
tion. 

We split up the K matrix in a diagonal Kq and non diagonal 
Xi] matrices, where parameter A is just introduced to account 
for the order in a perturbative expansion of the non-diagonal 
piece, 



We can then get the entropy Si{ m y . 

s l{m} * s l{m}A = f^(l- iog(o) 



k=l 



k=0 



(39) 



where the coefficients and tk are defined in Appendix A. A 
similar result for the single copy entanglement reads, 



Eic-J2 D) log (1 - + 0(l~ 10 ) 



1=0 



oo oo oo 5 

= E^)ET-£-E7S7 + °( rl0 )< 



k ^ ^ 11+3 
1=0 k=l 1=0 j=l 



(40) 



where Kj are the coefficients of the expansion given also 
in Appendix A. Finally, using Eq.( l34t and defining Tk = 

J2 J= o v 3 t k-3 and cr fc = J2j=o v 3 s k-j where v 3 are the co- 
efficients of the degeneracy expansion, we determine the con- 
tribution to the total entropy, for / = lo . . . oo, where Iq is big 
enough such that the approximations are valid , 

A ^-E^(c(6-^ +J )-E7g^Tj) 



$> c'(6-£+j)+E 



i=i 



16-D+j 



(41) 



where £(n) is the Riemann Zeta function and (,'{n) its deriva- 
tive. Defining Afc = v 3 K k-j, the single copy entangle- 
ment becomes, 



A^^E A ^fc(6-^+i)-E™ 

3=0 \ 1=1 



(42) 



The above results show that the sum over angular momenta 
I converges for D < 5. A radial discretization of a scalar field 
theory produces finite results for D < 5 and needs further 
regularization in orthogonal (angular) directions to the radius 
in higher dimensions. We will come back to this question 
later. 



K = K + Xri . (37) 

This expansion is somewhat tedious and non illuminating. 
Technical details are presented in Appendix A. The main ob- 
servation is that the first contribution i = 1 out of every set of 
£,i,{m},i elements is relevant and it can further be expanded as 
a series in Z _1 , 

e^,M,i = 7iE! + 0(r 10 ) (38) 



D. Area law scaling 

The analysis of the scaling law obeyed by the geometric 
entropy proceeds as follows. The analytical treatment of the 
chain of oscillators lead to the final sum over angular momenta 
in Ea.J33>. The computation of this sum requires polynomial, 
rather than exponential, effort as the size of the system in- 
creases. This justifies why large systems are accessible within 
this approach. The eigenmodes ^{ m j i are obtained by diago- 
nalization of matrices of order less than N. Finally the tail of 
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the sum over angular momenta is computed using the asymp- 
totic expressions given in Ea. J41l >. 

We have computed the geometrical entropy and the single- 
copy entanglement for different dimensionalities of the sys- 
tem. Within the range 1 < D < 5 we do observe the expected 
area law scaling 



S = k s (n,D,a,N) 



R 



D-l 



(43) 



as well as a similar scaling for the single copy entanglement 

D-l 



Ei = k E {»,D,a,N) ( - 



(44) 



where in all our considerations the lattice spacing can be taken 
a = 1. Fig [2 shows this perfect scaling for both measures of 
entanglement. 




100 200 300 400 500 600 700 800 900 1000 

Figure 1 : The entropy S and the single copy entanglement E± result- 
ing from tracing the ground state of a massless scalar field in three 
dimensions, over the degrees of freedom inside a sphere of radius R. 

The explicit pre-factor in the area law is regularization de- 
pendent but can be computed and compared with previous 
analysis. Fig|2 shows the result obtained for this pre-factor 
in the area law for the case of D = 3 and /i = as the size of 
the system increases. 

Good stability is already reached for N = 600, where we 
recover the result of |01 and complete it with the single copy 
entanglement 



= 0,D = 3,N 
k E (ji = 0,13 = 3, N 



oo) = 0.295(1) , 
oo) = 0.0488(1) 



(45a) 
(45b) 



Let us note that the ratio of the area law pre-factors for the 
entropy and the single copy entanglement is close to 6. This 
value is much larger than the factor of 2 computed to be the ex- 
act ratio in one-dimensional critical systems We thus con- 
clude that the amount of entanglement that can be extracted 
from a single copy of a system as compared to the asymptotic 
value for infinite copies does decrease with the dimensional- 
ity. 

We can analyze in more detail the dependence of our two 
measures of entanglement as a function of the dimensionality 



0.2963 | 

0.2962 - 

0.2961 - 

0.296 - 

0.2959 - 

0.2958 - 

0.2957 - 

0.2956 - 

0.2955 - 

0.2954 — 



0.0494 

0.0493 - 

0.0492 - 

0.0491 - 

0.049 - 

0.0489 - 

0.0488 — 



]l)l>2()()30l)4O()5O()6n(] 







300 

N 



Figure 2: Coefficient for entropy the area law in D = 3 as a function 
of size of the system. Good stability is reached for TV = 600. In 
the inset, the corresponding coefficient for the single copy area law 
is plotted. 



of the system. This is shown in Fig|3] for an TV = 60 and 
5 < n < 30 as a function of R ^ 1 and we verify that the area 
law is observed for any value of the dimension D. 

The robustness of area scaling law for arbitrary mass /j, is 
also readily checked (see Fig|3. The appearance of a mass 
term in the Hamiltonian produces exponential decays of cor- 
relators but does not affect the short-distance entanglement 
which is ultimately responsible for the area law. This sup- 
ports the idea that geometric entropy comes from the local 
neighborhood of the surface separating the region which is 
integrated out. The exponential decay of massive modes is 
immaterial and their contribution to the entanglement entropy 
is as important as the one coming from massless modes. 

Let us concentrate briefly in the dependence of ks and k E 
on the dimension D. Those coefficients present divergences at 
D = 1 and D = 5 (see Fig. [3]). The first one is due to the fact 
that in one dimension the strict power area law breaks down, 
since the limiting case carries a logarithmic dependence. For 
D > 5, as we have shown before, the sum over partial waves 
does not converge. This is due to the fact that we have regu- 
larized the Hamiltonian using a radial lattice. This regulariza- 
tion is insufficient to handle higher dimensional modes due to 
the increase of degrees of freedom per radial shell. To avoid 
this problem, a more elaborated regularization of the initial 
D-dimensional Hamiltonian is required. Such a regulariza- 
tion will likely have to break the rotational symmetry and will 
make the computations rather involved. 

We observe in Fig|3]that the entropy to single-copy entan- 
glement ratio verifies the expected limit 2, for D tending to 
1. 



E. Vacuum reordering 

Area law implies that entropy grows with the size of the 
system, that is, the eigenvalues of the density matrix, properly 
sorted from the largest to the smallest, decay in a slower way 
for larger systems. It has been numerically shown in Ref. 14811 
that this order relation between systems of different length 
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T 




I I'" —— I 4 1 1 1 1- — -" I 

1 1.5 2 2.5 3 3.5 4 4.5 5 

D 

Figure 3: Dependence of the geometric entropy and single copy en- 
tanglement slopes, ks and ks, on the dimension D for a massless 
scalar field. Note the divergence at D = 5 due to the insufficient 
radial regularization of the original field theory. 




1 1.5 2 2.5 3 3.5 4 4.5 5 

D 



Figure 4: Evolution of entropy to single-copy entanglement ratio 
S/Ei as a function of the dimension D. The line starts at a value 
of 2, as demonstrated analytically in 0] and grows monotonically. 
The higher the dimension is, the less entanglement is carried by a 
single copy of the system as compared to many copies. 



verifies the strong condition of majorization, a fact proven an- 
alytically for conformal field theories in Ref . l49ll . As the size 
of the system increases from L to L' > L, it is verified that 
Pu -< PL, where p^ and pi> are the set of eigenvalues for the 
corresponding reduced density matrices. 

Majorization relations characterize strong ordering. Every 
eigenvalue changes in a way that is consistent with a set of 
majorization constraints. We shall refer to this fact as vacuum 
reordering. 

We show that the same underlying reordering of the vac- 
uum is present in any number of dimensions. Unfortunately, a 
similar analytical treatment to the D = 1 case is out of reach 
because the conformal group in D > 1 is spanned by a fi- 
nite number of generators. As a consequence, there is no full 
control on the partition function of conformal field theory in 
D > 1 dimensions, which could be used to generalize the 
one-dimensional theorem. 

Vacuum reordering can be treated within our semi numer- 



ical approach. From Eq.( I28> we see that the reduced density 
matrix of the exterior of a ball of radius R, can be expressed 
as a tensor product of simpler density matrices, 

(N-n \ 
J| Pl{m},i{R) j , 
i=l / 

(46) 

where pi{ m } is what we call p out in Sec. dill Bl and pu m },i are 
defined in the same section. A similar composition applies for 
another size R', 

(N-n' \ 
n PIMA*?) ■ 
1=1 J 

(47) 

It is shown as a lemma in Ref.j48| that, if majorization rela- 
tions are satisfied by each pu m \ (R) and p;{ m j (R'), they will 
be also satisfied by p(R) and p(R'). Note, though, that it is 
not possible to follow the same argument for pu m \ i (R) and 
Pi{m},i(R') since n ^ n' . To make dimensions agree, we 
need to complete with identity operators the smallest set. We 
then find that some majorization relations for the subparts are 
obeyed in one sense, and the rest in the opposite one. Thus, we 
construct the density matrices p/{ m } (i?) and pi{ m } (R') doing 
the tensorial product of their components which are generated 
using Ea. J31> . Once we have their eigenvalues we are ready 
to check that if R < R', then 

Hout Pout (R) , (48) 

which means by definition 

k k 

I>'i<I> V* = l,...,oo (49) 

»=1 i=l 

where pi and p[ are the eigenvalues of p ou t{R) and p ut{R') 
respectively. For the I ~ N case, we have done a numerical 
computation with N = 60 and truncating the vector of eigen- 
values at the 50th element. Several dimensions D and traced 
sizes n have been studied, and all majorization relations are 
satisfied in all of them, as expected. When I 3> N, we can 
use the analytical results of the Appendix A to check the same 
result. 



IV. ENTANGLEMENT LOSS ALONG RG TRAJECTORIES 

We shall now exploit the control achieved on the eigenval- 
ues of the reduced density matrix in D dimensions to study 
how entanglement evolves along renormalization group trans- 
formations. This was studied for the quantum Ising model in 
Ref. 0. We shall now add equivalent results for the set of 
harmonic oscillators in D dimensions. Results will turn out to 
be qualitatively similar, reinforcing the concept of entangle- 
ment loss along RG flows. 

The renormalization of a bosonic field is particularly sim- 
ple since the Hamiltonian only carries one coupling, namely 
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the mass term. After a block transformation, the rescaling of 
fields is used to make the kinetic term to be normalized to |. 
The RG flow of the massive scalar field reduces to an effec- 
tive change of the mass. That is, the study of the long distance 
behavior of a correlator is viewed as taking a larger mass for 
the field, modulo a scaling factor. This implies the existence 
of two fixed points which are p = (ultraviolet, UV) and 
p = oo (infrared, IR). Since no other fixed point is possible, 
the RG flow must be monotonic in p. 

Entanglement loss comes along this flow. First, we study 
this change from a global perspective. We observe the obvious 
global loss of entanglement. For p = 0, geometric entropy 
grows with a slope ks(D, p = 0) for the massless field and it 
is zero for the p = oo case. Thus, 

Suv>S IR VR (50) 

This result is related to the c-theorem as discussed in Refs. 
SHU EH IHlliil], which states global irreversibility in the 
RG trajectory which interpolates between UV and IR fixed 
points. 

On top of this global loss of entanglement, the geometric 
entropy obeys a monotonic decrease along the RG flow. This 
behavior is illustrated for D = 3 in Fig|5]the entropy for dif- 
ferent masses where it is seen that 

p' > p k s (p') < ksifJ,)- (51) 

Thus, the system is more ordered as the mass increases. 
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Figure 5: Geometric entropy S for a sphere of radius R in D = 3 as 
a function of the mass p. Note that larger masses produce a smaller 
coefficient in the scaling are law. 

It is natural to pose the question if this order relation verifies 
also stricter majorization relations, that is, vacuum reordering. 
Specifically, we analyze whether p(p') and p(p), the density 
matrices corresponding to the free bosonic model with masses 
p' > p respectively, obey p(p') and p(p). 

p(p) -< ftp!) . (52) 

Using similar arguments as in the previous section, we only 
need to check that each pi{ 7n ] t i{p) majorizes pi{ m }, iif^')- 



Suv 

\Suv > Sir 
p -< p ' if p < p' 

Sir 

Figure 6: Entanglement loss along the RG trajectories seen in the 
space spanned by the eigenvalues of the reduced density matrix. 

Considering Eq.d3ll. that means, 

k k 

£(1 - Of < £(1 - Ot Vfc = 1, . . . , oo , (53) 

i=l i=l 

and therefore, 

(l-£ fe+1 ) < (l-£' fc+1 ) Vfc = l,...,oo. (54) 

This happens if and only if £' < £. As in the previous section, 
we have verified this fact numerically in the / ~ N regime, 
and analytically using the perturbation calculus done in Ap- 
pendix A. 

It should be noted that monotonic loss of entanglement is 
mandatory in such a simple model with a single parameter (p) 
controlling the flow. It is far from obvious that such entropy 
loss is rooted in a such a subtle reordering of the vacuum as 
the one dictated by majorization. 



V. CONCLUSIONS 

Area law scaling for the geometric entropy is present in har- 
monic networks of arbitrary dimensions. This follows from a 
computation that makes use of and analytical approach ca- 
pable of making an analytical extension of the computation 
to arbitrary D, followed by a final numerical resummation of 
angular momenta, whose tail is controlled analytically. 

A similar scaling law is observed for the single-copy en- 
tanglement. This result suggests that entanglement, whatever 
measure we use, scales with an area law due to the fact that en- 
tanglement is concentrated on the surface of the region which 
is traced out. The ratio of single-copy entanglement to geo- 
metric entropy tends to zero as the dimension of the network 
increases. 

It is natural to interprete a change in the size of the sub- 
system which is traced out as well as any modification of 
the parameters in the Hamiltonian as a probe on the vacuum. 
Our explicit computations unveil ubiquous vacuum reorder- 
ing governed by majorization relations of the vacuum state re- 
duced density matrix eigenvalues. Geometric entropy scaling 
is just one manifestation of this set of order relations. 

The fact that finite PEPS support an area law scaling makes 
them a natural tool to investigate regularized quantum field 
theories. 
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Appendix A: PERTURBATION THEORY 

We need to perform a perturbative computation for large 
momenta in order to determine the contribution to the total 
entropy and single copy entanglement of all angular momen- 
tum modes. 

We organize our computation in three parts. In the first part, 
we carry out perturbation theory with matrices, following the 
same steps as explained in sec. JIII B I when considering the 
aproximation I ^> N. This will produce an analitical expres- 
sion for the £'s parameters. The second part of the computa- 
tion consists in Taylor expanding the above results for £ in a 
series in Finally, we will get the entropy and single copy 
entanglement contributions, expanding the entropy and single 
copy modes in terms of l^ 1 powers, and summing over I. In 
this sum I take values from Iq until infinity, where Iq must be 
sufficiently large, such that all aproximations done previously 
are right. 



1. Computation of the £ parameter 

Let us recall that, for I S> N, the non diagonal elements of 
K in Eq. (I23> are much smaller than the diagonal ones. That 
gives us the possibility of setting up a perturbative computa- 
tion. 

We split up the K matrix in a diagonal Kq and non diago- 
nal A77, matrices where A is just introduced to account for the 
order in a perturbative expansion of the non-diagonal piece, 



K = Kr 



Ar? 



(Al) 



We will follow the steps described in Sec HUB I We expand 
fl = \[K in its different contributions to order A, 



n = fi + Ae + X 2 e + A 3 e + 0(A 4 ) . 
To get each term we impose the condition SI 2 = K, 

( e )ij = e i^i+l ,j + € j$j+l,i 



(A2) 



ef + e 



H+2,j 



'I I J 

(ee + ee) l 



where ilj and ej are defined since 



(A3) 



ft; 



U(l + D-2) 



1 + v 

3 + 



D-l 



D-l 



(A4) 



We structure fl in three matrices A, B and C, 

A = ft(l-=-n, l-=-n) = A a + XA X + X 2 A 2 + <3(A 3 ) 
B = ft(l-=-n,n + 1-hiV) = 

= XB + X 2 B 1 + X 3 B 2 + 0(A 4 ) 
C = ft(n + 1-t-N, n + 1+N) = 

= C + Xd + X 2 C 2 + 0(A 3 ) . (A5) 

From these matrices, we define /3 and 7 which we write in 
series of A 

/3 = ^B T A~ X B = A 2 /? + A 3 /3i + A 4 /3 2 + 0(A 4 ) 
7 =C - /? = ft + Ae + O(A 2 ) , (A6) 
where, 

1 £ 2 

A) = 2 B ° A o B o = 2jY~*i,i^',i 

p! = ^(BfA^B + Bf ) Av 1 B 1 + B^A^Bo) 

£ 2 e ^ +] 

2f2 ra f2 ra + f2„+2 ' 3 ' 3 
= \{BlA^B + BlA^Bz + B^A^B a 
+ BlAo 1 B 1 +B^A^ 1 B + B^A^ 1 B 1 ) . (Al) 

We shall see later, that at 2n order perturbation in A, only 
(/32)ii an d (^2)22 of (3% are necessary. Then, 



(#2)11 = 



n J n— 1 . 71— 1 > ?i 



ft. 



c n-l 



ze n-l 



ftn-1 (ftn+1 + ftri-l) 2 ftn-l(ft«+l + ft«-l) 



e n+l e n+l + e n 



ftn(ftra + ftn+1 ) \ ftn + ftn+2 2ft n _|_i 

,2 , 2 2 
n ^ c n— 1 1 n— 1 



257„ fin+1 + ^n— 1 



(/3 2 ) 



22 



c n c n+l 



2f2 n (f2 n + Qn+2) 2 
Let us diagonalize 7, 

1D = VjV T 



(A8) 



(A9) 



where V is an orthogonal matrix (VV T — 1). Therefore, the 
eigenvalues are 



N — 7 



dct( 7 -wl)= Y[ ( n n +l -w) + 0(X 2 ) = 

1=1 

= + 0(A 2 ) , 



and 



(7i?k = (fin +I + 0(A 2 ))r5, 



(A10) 



(All) 
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If we impose <|A9) over V = V + XV 1 + X 2 V 2 + 0(A 3 ), we 
obtain, 



V = l 



(Vi) 

(V 2 )n = 



13 ~ n 
1 



n+j 
e n+l 



Si+lJ + 



-n+3 



--rt 



(A12) 



Once we have V and -jd we are able to compute /?' = A 2 (/% + 
A/3J + A 2 /^ + 0(A 3 )), which is defined by, 

= 7d *W T 7 ~' . (A13) 

Thus, 

Pi = (loho [Pi + V (3i + fcV? + Vi0o + PoV?] (^)o 
& = (7 D ') 2 /3o(7^)o + (n D h )M D h )2 

+ (l D ho [P2 + Vifa + fa.V? + V 2 f3o + PoV 2 T 

+V 1 (3 V 1 T ] (7^)0 
and therefore, 



(A14) 



and c„ and e n are respectively (/3 2 )n and (0' 2 ) 22 . We can 
observe now, that if we had not found the second order contri- 
bution of (P')n and {fi') 22 , we would not have been able to 
compute the eigenvalues of /3' to this order. 
Diagonalizing j3' , we find the eigenvalues v\ and v 2 , 

v,=X 2 (a n + X 2 ( Cn + -)+0(X 3 ) 



« 2 = A 4 (e„-^l + 0(X 5 ) = + 0{X 5 ) , (A18) 



which allows us to compute the £j's parameters, 



and which read 



6 = y( a » + A2 ( c » + ^) 

6 = + O(A 5 ) 
6 = 0(A 7 ) Vi>2. 



2. Expansion of £ in terms of Z 1 powers 



(A19) 



(A20) 



W'oh 



-<5i,i<5j, 



(#) 



+1 



1 



£n+l 



X 



+ 



1 



£l n +l — Qn+2 + &n+2 



(<5i,l^,2 + <5i,2<5j,l) 
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(A) 
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c n+l 



,2 / 0,2 

ze n+l 



2^n^n+l \ (^ri+2 + ^ri)(^n+l _ ^71+2) 
1 ^ 2 "n+1 + ^n+2 2 ^« + 



2^22 



oo2 I "+ 1 
f 2 F 2 

t n t n+l I 



Qn+1 — Qn+2 



1 



2^n+2^?i \ — ^ri+2 "71 + " 



n+2 



(A15) 



It will be useful to write j3' in its matrix form, 
/ a n + X 2 c n \(L . . . \ 



/3' = A 2 



where, 



Xd n 




A 2 e„ 




+ 0(A 5 ) , (A16) 



2^71+1^?! 



1 



1 



^n+2 V — ^n+2 + ^ 



n+2 



(A17) 



We rename £1 as £, and neglect the rest since at this order 
they are and no contribute neither to the entropy nor to the 
single copy entanglement. We are interested in expanding £ 
in powers of l^ 1 . To do this, we have to expand first ftj and 



9 o« 
i=0 



o(r 9 ) 



(A21) 



where, 

n 

0< 2 

n 

n (4 

fil 5 



71 

D-2 

2n 



(D-2) 2 



2 8n 
(D - 2) 3 - 4(£> - 2)n 2 w„ 
16n 

5(,D - 2) 4 - 24(D - 2) 2 n 2 w„ + 16n 4 w 2 
128^ 

7(£) - 2) 5 
256n 

-40(£) - 2) 2 n 2 w„ + 48n 4 w 2 
25fe 



(A22) 



No more coefficients have been presented here since they have 
huge expressions and they don't shed any light on our argu- 
ments. Using Q n we can obtain the expansion of e„, 



1 6 

i=0 



(A23) 
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where, 



where, 



(</) 



(0) 



c (0) 



f ( 2 ) = f (°) 



n,n+l 



+ ft (0) 






"n+l 




"ri+l 






IU 0) 


+ r> (0) 

+ "ri+l 



so = Co - Co log Co 
Si = -Ci log Co 

d 2 



(2) , o( 2 ) 



S2 



S3 



2Co 



C2logCo 



& - 6C0C1C2 
6C0 2 



-C3log'Co 



n+1 



o (0) + n {0) 

tl \ J n+1 



(A24) 



and 



Once we have f2„ and e„ in series of I 1 , we can expand C, 

e = ^Ef + or 10 )> 

i=0 



t, = 4& < i < 3 

U = 4(C 2 + C 4 ) 

t 5 =4(2CoCi +Cs) • 



(A29) 



(A25) 



To determine the contribution to the entropy of all modes with 
the same I, we use the expansion of the degeneration, 



with 



v(l,D) 



l + D - 1 
I 



l + D -3 
I - 2 



(e£ 0) ) 2 



^" 2 E 



(A30) 



fc=0 



Ci = e 



(o). 



2e} : 



^fi^fi^ - ei 0) ff^^i+i + fi^fi^) which allows us to sum over all the possible values of {m}, 



(rf) 2 + (n%y 



(A26) 



{m} 



i=0 



Although C depends on the number of oscillators which we 
trace out, we have omitted the subindex n to simplify the no- 
tation. 



3. The entropy 

The contribution to the entropy of a (I, {m})-mode be- 
comes, 



(A31) 

where r k = V* t| /<,//. , and <j k = V, Vj8k ,. Finally, 
we can compute the contribution to total entropy, for I from 
Iq to 00, where Iq is big enough such that these aproximations 
are justified. 



3 \ i=i 



(A32) 



C; 



Si{ m } = - log (1 - CO - y~t^ logCz 

-E^-mc)V. 



fe=i 



(A27) 



If we substitute C, 



being the Riemann Zeta function, and C'(^) i ts deriva- 
tive. 



4. The single-copy entanglement 

We can do the same as for the entropy to find the contri- 
bution to the single-copy entanglement for large values of I. 
First, we expand the contribution to the total single-copy en- 
tanglement of the {I, {m})modes, 



1 ^ s k +t k \ogl n(1 - w . 



I* 



k=0 



I* 



(A28) (E 1 ) l{m} ~ - log (1 - CO = J2 jS? + °( rl °) . (A33) 



i=0 
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where, 

Ki = £i < i < 3 

«4 = (f + &) 

«6 = + (A34) 
Next, we sum for all possible values of {to}, using Eq. (IA30t . 

5 ^ 

{E 1 ) l = v{l,D){E 1 ) l{m} =Y,w^+ i , (A35) 

i=0 

where Afc = X^=o v i K k-j- Proceeding as before, we finally 
get 

El cJ2 Aj ( C(6 - + i) - £ ^zi+j J • (A36) 

7=0 \ 1=1 / 
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